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We derive a kinetic theory for a spin-1/2 Bose-condensed gas of two- level atoms at finite tem- 
peratures. The condensate dynamics is described by a generalized Gross-Pitaevskii equation for 
the two-component spinor order parameter, which includes the interaction with the uncondensed 
fraction. The noncondensate atoms are described by a quantum kinetic equation, which is a gen- 
eralization of the spin kinetic equation for spin-polarized quantum gases to include couplings to 
the condensate degree of freedom. The kinetic equation is used to derive hydrodynamic equations 
for the noncondensate spin density. The condensate and noncondensate spins are coupled directly 
through the exchange mean field. Collisions between the condensate and noncondensate atoms give 
rise to an additional contribution to the spin diffusion relaxation rate. In addition, they give rise 
to mutual relaxation of the condensate and noncondensate due to lack of local equilibrium between 
the two components. 



I. INTRODUCTION 



In this paper we present a kinetic theory for a trapped spin-1/2 Bose-condensed gas describing nonequilibrium 
collective dynamics of the density and spin at finite temperatures. Our work has many points of contact with several 
strands of research conducted over the past few decades. In a way, this article provides a conceptual bridge between 
the very active field of Bose-Einstein condensation (BEG) in dilute gases and earlier work done on spin-polarized gases 
in the 1980's PjQiIIISIEQj in which the field of ultracold gases has its roots. It is interesting that, while cooling down 
to BEG was a central goal of the sp in-p olarized hydrogen work, spiri waves in a nondegenerate gas became a major 
topic for both experiment [l|llp.[iallll and theory 13, lllIEIlEEHis, 19, 20]. After BEG was eventually observed 
in 1995 in the alkali atoms |2l]. 122. l23l| . however, the physics of spin waves in the new breed of ex per iments had been 
largely overlooked — until very recently, when spin waves were observed at JILA [3, HE HE HE HE HE HE Isills^. 
So an obvious direction for theory is to extend the earlier kinetic theories that were developed for a dilute spin-1/2 
gas above Tc 

[IEIIllIEIiEIIEIiEIIEH3 

into the Bose-condensed regime, which is the main thrust of our paper. 
Alternatively, our work can be viewed as an extension of the recent work on finite temperature kinetic theory for 
a single-component Bose-condensed gas [sE 13 HE H^ to now take into account the spin degree of freedom of the 
atoms. 

Before placing our work in a broader context, it is useful to first clarify the type of system we have in mind for our 
kinetic theory. Experiments at JILA [2J, HE HE HE HE HE HE HEl have explored various properties of a dilute Bose 
gas of two-level atoms confined in a magnetic trap at finite temperatures. Using a two-photon coupling field, ®^Rb 
atoms can be prepared in a superposition of the two hyperfine states [F — 1,Mf — —1) and (2, 1). This two-level 
atom can be treated as a spin-1/2 system. These states are particularly nice to work with because they have nearly 
the same magnetic moment, and thus experience the same trapping potential (to first order in the magnetic field 
gradient). In ^^Rb, spin-exchange losses in the magnetic trap are reduced by more than a factor of 1000 compared 
to other alkalis due to a fortuitous near-degeneracy of the singlet and triplet scattering phase shifts [43, HE HE HE| ■ 
Ramsey fringe spectroscopy of the hyperfine splitting shows that the internal coherence-or transverse spin-can be 
preserved throughout the gas for over a second, which is much longer than the coUisional relaxation time 0,H3- 
In this scenario, the transverse spin polarization can be treated approximately as a conserved quantity, although for 
long enough times the spins will dephase due to inhomogeneities. These properties make this system ideal to study 
as a model spin-1/2 Bose gas. 



A. Pre-BEC work on spin waves 

Gollective spin behavior is a well known property of metals with ferromagnetic or antiferromagnetic ordering |48j . 
The idea that spin waves could also propagate in a paramagnetic system, such as '^He or a normal metal, was first 
pursued in 1957 by Silin ^EHEIj who generalized the Landau Fermi liquid theory to treat the effect of an external 
magnetic field. By including the off-diagonal magnetization terms in the single-particle distribution function, Silin 
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derived a coUisionless spin kinetic equation and showed that transverse spin waves could propagate in a finite magnetic 
field. The first experimental verification of spin waves in a paramagnetic system was reported by Schultz and Dunifer 
in 1967 by measuring the electron spin resonance spectrum of thin metallic slabs Shortly thereafter, Leggett 

and Rice predicted that, as a consequence of the collective spin behavior, the spin diffusion coefficient should have 
a peak as a function of temperature, which could be measured using an NMR spin echo technique (53i. .54] . This 
so-called "Leggett-Rice effect" was then verified in 1972 by Corruccini et al. |H3 in '^He-^'He mixtures and normal 
liquid ^He. A more direct observation of spin waves in liquid "^He came much later when the frequency spectrum of 
standing spin wave modes was measured by Masuhara et al. in 1984 ,56, 57]. 

All of this work dealt with dense Fermi liquids, so it came somewhat as a surprise when two independent predictions 
appeared in 1982 by Bashkin and Lhuillier and Laloe 0, 0] that transverse spin waves can also propagate in 
dilute nondegenerate Bose and Fermi gases. A few years earlier, dilute atomic hydrogen gas had been stabilized 
against molecular recombination at a few hundred degrees millikelvin by polarizing the electronic spin in a high 
magnetic field [HI, [s^ |63, which provided a perfect testing ground for these predictions. In order for a dilute gas 
to exhibit collective spin behavior, the temperature has to be low enough that the thermal deBroglie wavelength Ath 
is larger than the effective range of interaction ro, i.e., Ath > ''o- If this condition is satisfied, then identical-particle 
symmetrization of the scattering wavefunction for two colliding atoms is required and this gives rise to the quantum 
exchange interaction needed for collective spin behavior. The dilute hydrogen gas certainly satisfied this criterion, 
but yet was still far from the quantum degenerate regime since n^/'^Ath ^ 1, which means that quantum statistics do 
not play a role in thermodynamic properties. To help dispel this paradox, Lhuillier and Laloe offered the following 
picture, which emerged from their microscopic treatment [liL Il5|: When two colliding atoms scatter off each other, 
their individual spins precess about their net spin due to identical-particle symmetrization, even if the interaction has 
no explicit spin dependence. The cumulative effect of these spin rotations in successive collisions throughout the gas 
is a collective oscillation of the transverse spin density. This "identical spin rotation effect" manifests itself in the 
kinetic theory as a nonlinear mean field term proportional to the local spin density S(r,t), which exerts a torque on 
the spin current J(r, t). The sign of this term for a Fermi gas is opposite to that of a Bose gas. 

The first experimental verification of collective spin oscillations in a dilute gas came in 1984 by Johnson et al., who 
measured the spin wave resonances in the NMR spectrum of dilute spin polarized hydrogen gas in the hydrodynamic 
regime y.[lfij . Later on this work was extended by Bigclow et al. to study the system in the low-density coUisionless 
regime [ll|. Shortly after the first hydrogen experiment, similar studies were done by Nacher et al. with dilute spin 
polarized '^He gas [a0] and by Gully and MuUin with a solution of ^He in "^He, which can be described as a dilute 
nondegenerate gas of quasiparticles. 

Later work on collective spin effects in paramagnetic systems focused on understanding the effects of quantum 
degeneracy. In 1988, Jeon and MuUin |12, 18, 2Q] and, independently, Ruckenstein and Levy jf^, used the Kadanoff- 
Baym nonequilibrium Green's function approach [Glj to obtain the kinetic equations describing quantum-degenerate 
(but nonsuperfluid) spin- 1/2 Bose and Fermi gases. The focus of these studies was to confirm with a more rigorous 
theory an earlier prediction made by Meyerovich in 1985 JJ^ that in a degenerate Fermi system the spin-diffusion 
relaxation time t± for transverse spin waves should differ from the longitudinal relaxation time r|| . Jeon and MuUin 
showed that, while the longitudinal relaxation time exhibits the expected divergent behavior T|| ~ l/T'^, the transverse 
relaxation time saturates to a finite value as T — > 0. In addition to dependin g on temperature, this anisotropy also 
depends on the degree of polarization (which is why the earlier theories of Silin |49ll50| and Leggett overlooked 
it, since they both take the limit of zero polarization in equilibrium). Experimental evidence supporting this prediction 
appeared in the mid 1990's from Wei et al. and Ager et al. |63, who used NMR spin echo techniques to measure 
the spin-diffusion coefficients in spin polarized ^He. Very recently, however, similar experiments were conducted by 
Vermeulen and Roni |64j at a much higher polarization and they find no evidence for polarization induced spin wave 
damping |^. It is clear that further experiments at higher polarizations and lower temperatures could give more 
definitive results. Two-component degenerate Fermi gases of alkali atoms, which are being actively studied by several 
different groups j93], could play a deciding role in answering this question, since polarizations of 100% are possible in 
those systems. 

The effects of degeneracy in a spin-1/2 Bose gas were not explored in great detail before BEG was achieved in 
1995, since ex peri ments had not operated in this regime before that time. Sparked by the paper in 1976 by Stwalley 
and Nosanow 5,8] that made a strong case for BEG in spin-polarized hydrogen, a number of theory papers appeared 
over the following several years on BEG in a dilute spin-1/2 gas (see for example, Ho |6^ and references therein), 
emphasizing that the internal degrees of freedom made the system potentially much richer than superfluid *He, and 
would have a closer resemblance to superfluid '^He in terms of its superfluid dynamics and symmetry breaking 
properties|i62j. The theory of spin waves in a Bose-condensed spin-1/2 gas, however, was not pursued, despite the 
fact that predictions for a classical gas by Bashkin ^3] and Lhuillier and Laloe 0, llSB app eared around that time in 
the early 1980's. Steps were certainly taken in this direction when Jeon and MuUin^j^i 13 113 and Ruckenstein and 
Levy (19j] extended the kinetic theory to lower temperatures where quantum statistics play a role , but spin waves 
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in the Bose-condensed regime were so far never treated. 

B. Post-BEC work on spin waves 

After BEC was achieved in 1995 with magneticaUy trapped alkah atoms, there was no immediate attempt to 
investigate spin waves in the Bose-condensed regime. In general, it is very difficult to maintain a spinor gas in 
a magnetic trap, thus the main focus of research was on single-component condensation, in which only a single 
hyperfine state is occupied. Unlike in the earlier hydrogen experiments, the alkali atoms are typically placed in a very 
weak magnetic field (around a millitesla), and so the total hyperfine spin governs the magnetic confining potential. 
There are only three low-field seeking states in ^^Rb, ^^Na, and ^Li: {F = 2, Mp = 2), (2, 1), and (1, —1), which can 
be magnetically trapped. Conventional wisdom is that having a true spinor system in a magnetic trap is precluded by 
spin-exchange collisions, which transfer atoms to untrapped hyperfine states. Typical loss rates due to this density- 
dependent process are on the order of 10~^^cm'^/s for the alkali atoms, which can limit the lifetime of multicomponent 
gases to less than a millisecond. One way to circumvent this problem is to load the atoms into an all-optical trap, as 
first done at MIT in 1998 with ^^Na atoms prepared in the F = 1 hyperfine multiplet f69„.70.j . which are simultaneously 
trapped in the same confining optical potential. In this spin-1 system, spin-exchange collisions play a beneficial role by 
allowing the relative populations between the spin states to equilibrate in order to reach the thermodynamic ground 
state. 

To have a spin-1/2 Bose gas, we need two hyperfine states of atoms in a gas. In Appendix IXI we elaborate on 
several two-level systems at our disposal in dilute atomic gases. As mentioned above, it turns out that in ^^Rb, 
because the singlet and triplet scattering phase shifts are nearly degenerate, the spin-exchange loss rate is on the 
order of 10~^^cm^/s, so multicomponent ^^Rb gases can have lifetimes of more than a second |^ |3 EE By 
preparing ^^Rb atoms in a superposition of the hyperfine states (2,1) and (1,-1), this two-level system in many 
respects resembles spin polarized hydrogen, in which the two states correspond to the up and down states of the 
nuclear spin. Conceptually, this system is a cross between a spinor gas and a binary mixture of distinct atomic 
species. Under many circumstances, it can be appropriately treated as a binary mixture because, in the absence of an 
external coupling field, there are no intrinsic mechanisms that lead to interconversion of the two spin states, so that 
the relative spin population is conserved Unlike in a mixture of two distinct atomic species, however, the relative 
phase between hyperfine states can play an important role in the collective dynamics due to the exchange interaction 
discussed above. At T = 0, the exchange interaction is absent from the coupled Gross-Pitaevskii (CP) equations 
describing the two-component condensate, and so the system appears equivalent to a binary mixture 95]. At finite 
temperatures, however, due to the exchange mean field, the relative phase between the components can strongly affect 
the collective dynamics of the relative density, which makes the system absolutely distinct from a binary mixture. Put 
another way, a spin-1/2 system has an extra degree of freedom — the transverse spin — which is absent from a binary 
mixture. 

This point may not have been full y appr eciated before the recent observation of spin waves at JILA [23. |2^ 
because all previous experiments [s^, IsM l39l liof and most theoretical treatments of the two-component ^^Rb gas 
had concentrated mainly on zero temperature behavior of the condensate. In contrast to these earlier works, the 
recent JILA experiments were done in the high-temperature regime T 2Tc, where the quantum degeneracy has 
little effect on the thermodynamic properties of a gas. All the atoms are initially prepared in the same (1,-1) state 
(i.e., all in the spin- up state), and then a 7r/2 pulse is applied to tip the spins into the transverse direction. The spin 
vector then precesses about the longitudinal axis at a rate proportional to the energy difference between hyperfine 
states. Due to the mean field and differential Zeeman effects, the local frequency splitting between hyperfine states 
varies approximately quadratically with position. This inhomogeneity initiates collective spin dynamics through the 
exchange mean field. A large inhomogeneity induces nonlinear spin oscillations, which cause the striking spin-state 
segregation initially observed p^ . This inhomogeneous frequency splitting can be also made arbitrarily small to 
study the linear response of the system. This technique was used to probe intrinsic collective spin oscillations |28j |. 
The collective spin dynamics observed at JILA above Tc are well understood by the theory based on the Boltzmann 
kinetic equation with spin degree of freedom |2^ |2^ [SJ, ■ 

One exciting advancement of the recent spin wave experiments at JILA '24*, "2^ is the technique to obtain spatially 
resolved images of spin dynamics in a gas. The density profile of either spin state is measured using absorption 
imaging. Together with the Ramsey fringe technique, integrated spatial profiles of the longitudinal and transverse 
phase angles of the local spin density can be extracted from experimental data, as shown in the stunning images in 
McGuirk et al. [2^ . This is in sharp contrast to the earlier experiments on the spin-polarized gases, where NMR was 
done on the system to obtain the frequency spectrum spatially averaged over the entire sample. This new technique 
stimulates theoretical investigation of the spatial structure of the spin dynamics p9l | that can be directly compared 
with experimental observations [28|. 
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C. Outline of this paper 



In this paper we derive coupled dynamical equations for the condensate and noncondensate components, both of 
which have spin-1/2 internal degrees of freedom. For a single-component Bose-condensed gas, finite-temperature 
kinetic theory has been a rnajor field of theoretical study, and several groups have derived kinetic theories using 
different approaches |3a. [SJ, [S^, ■ Among them, our present theory has the most overlap with the kinetic theory 
by Zaremba, Nikuni, and Griffin (ZNG) '3^, which uses the semiclassical description of the noncondensate atoms 
within the Hartree-Fock approximation — a theory closely related to the pioneering work of Kirkpatrick and Dorfman 
[tH- Recent numerical simulations by Jackson and Zaremba [t^ showed that the ZNG kinetic theory provides an 
excellent quantitative description of various dynamical phenomena in trapped Bose gases at finite temperatures. It is 
thus natural to also expect the same kind of semiclassical kinetic theory to be effective for a spin-1/2 Bose-condensed 
gas. In Sectional we give a detailed derivation of a generalized GP equation for the condensate and a semiclassical 
Boltzmann kinetic equation for noncondensate atoms. These two components are coupled through the Hartree-Fock 
mean field as well as collisions (corresponding to "C12" collisions in the ZNG theory). In the ZNG theory, C12 
collisions play an important role in coUisional damping of the condensate collective modes [Tsj, condensate growth 
[73. Its! , and vortex nucleation [t^ . We will show in this paper that the condensate-noncondensate collisions also play 
an important role in the spin-1/2 gas in bringing the two spin components into local equilibrium with each other. 

We emphasize that the dynamics of the spin-1/2 Bose-condensed gas at finite temperature is potentially more 
interesting than that of a single-component system because the transverse spin of the thermal cloud can have strong 
collective dynamics even in the very dilute limit. In contrast, the mean field of the noncondensate plays a very minor 
role on the density fluctuations in a single-component Bose gas, (mainly as a source of damping of the condensate 
excitations) |35, 72]. To gain more physical insight into the coupled spin dynamics, in Section lTTll we recast the kinetic 
theory derived in Section ^] into a form of spin hydrodynamic equations, which are written in terms of the spin and 
spin current for the condensate and noncondensate components. At finite temperatures, the spins of the condensate 
and thermal gas interact strongly through the exchange mean field. The role of the exchange interaction between the 
two components has been discussed by Oktel and Levitov in a spatially homogeneous system. Our hydrodynamic 
equations also involve the collisional relaxation term due to "Ci2-type" collisions, which try to bring the thermal 
cloud and condensate spins into local equilibrium. 

In Section IIVI we investigate the effect of different scattering lengths on the relaxation of the transverse spin 
polarization. A similar discussion was given in a recent paper by Bradley and Gardiner |3lj | . In Section ^ we 
summarize this paper, and discuss some possible future applications of our kinetic theory. All the detailed calculations 
of the collision integrals and transport relaxation times are given in the Appendices. 



II. DERIVATION OF THE SPIN-1/2 KINETIC EQUATIONS 



We consider a trapped Bose-condensed gas of atoms with two internal states, which are denoted as |1) and |2). For 
example, these can be any of the systems shown in Fig. 1 as discussed in the Appendix. The state of a single atom is 
given by a spinor wavefunction 



= (^i(M)) ='^l('■'*)|^>+'^2(^'^)|2^■ 



Bold-faced variables are vectors in coordinate space and a single underline indicates a spinor variable. The single-atom 
Hamiltonian is given by 



fio(r) 



2m 



V2 + C/e,t(r) 



l + |f^(r)-£, 



(2) 



where the double underline indicates a 2 x 2 matrix, and the vector symbol (e.g. fl) indicates an SU(2) spin vector. 
The first term describes the center-of-mass motion of an atom in a harmonic trap 



TT /\ ^ (22, 22, 2 2\ 

C/ext(r) = -m(w^a; -t- w^^?/ +u^z ) 



(3) 



The second term describes the evolution of the spin state of the atom. The dot product can be expanded as 

l1(r) • £ = f7^(r)£^ -f ^^(y)^ -f ^^^(r)£^ (4) 



5 



F=2 



Magnetic trap: ^^Rb 

(F=2, Mp=2) 

f(2,^)\ 

(2,0)' 



(2,-2) 



F=1 




two-photon 
transition 



Optical trap: Ti, ^^Na, ^^Rb 




F=2 



B 



FIG. 1: Several possible two-level systems. Panel A shows the two states used in the JILA experiments, where Rb atoms 
are held in a magnetic trap. Panel B shows two different cases that would make good two-level systems for the alkali atoms 
confined in an optical trap. 



The matrices {j = x, y, z) are the usual Pauli spin matrices. The quantity 0(r) describes an external field (analogous 

to an external magnetic field in a true spin- 1/2 system). The transverse components J7^(r) and 17^ (r) describe the 
possible coupling between the two spin states (for example, if a microwave field or laser is applied that is tuned close 
to resonance) and the longitudinal component VL^(y) describes the energy offset between the two internal states. In 
addition to being position dependent, the harmonic trap f/cxt and the external field 51 may also vary in time. We 
describe the collisions between two atoms in the gas by the contact interaction term 



V^2{r,v') = (5(r-r') Vb 1 1 + (1 £^ + £^ (g) 1) + Vz;^£^ 



(5) 



where (g) is the tensor product. The subscript ®2 denotes a 4 x 4 matrix formed from the tensor product of spin 
matrices and Vq = {gn + 522 + 2gi2)/4, = (.gn - .g22)/4, and V^z = (511 + .922 - 25i2)/4. The tensor V^2 has the 
following four components: 



{i\{j\V^2m)^grj, (»,j = l,2). 
The interaction strengths for the various collision processes have the form 



9ij 



(6) 



(7) 



where a^- are the s-wave scattering lengths describing collisions (e.g. \i)\j) The scattering lengths for a 

specific two-level system must be obtained from a detailed multichannel scattering calculation. For the case of *^Rb 
prepared in the states |1) = \F = 1,Mf = —1) and |2) = \F = 2,Mf = 1), these scattering lengths have been 
determined to be an = 100. 9ao, 012 = 98.2ao, a22 = 95. Gap, where oq is the Bohr radius 1^. The near degeneracy of 
these scattering lengths gives the result Vzz ~ in Eq. © and the first term proportional to Vb dominates over the 
second term. If one makes the approximation of equal scattering lengths, then there is no explicit spin dependence in 
the interaction Vz — Vzz — 0. 

In quantum field theory, the many-body Hamiltonian of this system is given by 



E 



dr^J(r)(z|^o(r)|j)V^,(r)-H^^ 



dripj (r)-0j (r)V'j (r)'0, (r) , 



where ipiir) is a Bose field operator satisfying the commutation relation 



(8) 



(9) 



The hat indicates a second quantized operator. An important property of this Hamiltonian is that, in the absence of 
the transverse coupling field (i.e. il^ = il^ = 0), the total number operator for each state, given by 

(10) 
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commutes with H, which means that the population of each state is separately conserved. This property, which also 
occurs in a binary mixture of distinct atomic species, means that the two levels will not become thermally populated 
in equilibrium (in our model, there are no state-changing collisions that lead to a net interconversion between states). 

We next define the time evolution of the system in terms of the statistical density operator, and introduce the 
Heisenberg representation for the field operators. The dynamics of the system is described by the density operator p{t), 
from which one can obtain the expectation value of an arbitrary operator O (which has no explicit time dependence) 

{d)t = trp(t)d. (11) 
The state of the many body system evolves in time according to 

^h^^p{t)^[H,m]■ (12) 
With Eq. (|ll|l and Eq. H12|l . the equation of motion for the quantity {0)t is given by 

inj^{d), = {[d,H]),. (13) 

It is convenient to introduce the time evolution operator U{t, to), which obeys the equation of motion 

ih-^U(t,to) = HU(t,to), (14) 
dt 

with U(ta,t()) = 1. Here Iq is the time at which the initial nonequilibrium density matrix p{to) is specified. One can 
then express the time evolution of p{t) as 

p{t)^y^{t,to)p{to)UHt,to)- (15) 

Thus, the time evolution of {0)t can be written as 

{d)t = tTU{t,to)p{to)UHt.to)0 = tTp{to)U\t,to)dUit,to) = (6(t)), (16) 

where 0{t) = W {t,tf))OlL{t,t()) is the operator in the Heisenberg picture, which obeys the Heisenberg equation of 
motion 

^n^ = [dit),H{t)]. (17) 

The Heisenberg equation of motion for the Bose field operator is then given by 



at 



2m 



+ C/cxt(r) 



V3.(r, t)+J2 .9'u^](i-, t)iPj{r, t)Mr, t) 



-^|(^|^1(r,^).£|J)^,(r,^). (18) 



The key idea in capturing the physics the Bose-Einstein condensation is to introduce the broken-symmetry order 
parameter 

t(r,t)= f^t;)) =<i>i(r,i)|l)-Fa>2(r,t)|2), (19) 



,$2(r,t) 
with 

$,(r,i) = (^,(r))t = (V;,(r,i)). (20) 

In addition to the usual U(l) gauge symmetry breaking, the order parameter defined by Eq. (|19|l also breaks the 
SU(2) symmetry, which is associated with the spin-1/2 internal degree of freedom. The noncondensate field operator 
is defined by 

^A»(r,i) = V',(r,t)-$,(r,t), (21) 
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with {'ip(r,t)) = 0. This noncondensate operator satisfies the equal-time commutation relation 



The equation of motion for $i can be obtained by taking the expectation value of Eq. (|18|l . which yields 



(22) 



El 



2y72 



2m 



+ C/c(r,t) 



<5y+ 2(«|f^c-£|j)<{ $j(r,i) 



[my(r,t)$*(r,t) + (^](r, i)?^,(r, t)i^,(r, i)) 



(23) 



Here we define the effective potential Uc and coupling field f2c, which include the mean- field interaction of the 
condensate with itself and with the noncondensate 



Uc = C/cxt + ^[511"! + 522"-2 + gi2{ni + 712) + giihi + 322^12] 



(24) 



.912 ^ - 

— — Lian, 

h — 



where "Tr" denotes the spin trace. The condensate and noncondensate densities are defined by 

(«|Bp(r,i)|j> = ney(r,i) = $*(r,i)$i(r,i), 



(25) 



(26) 



(«li(r,i)|.?) = n^J{v,t) = (V;j(r,i)V'.(r,i)) 



(27) 



with rici = ricii and hi = ha. The total density of internal state i is = rici + hi. The second term in Eq. H25|l directly 
couples the spins of the condensate and noncondensate and is present even when the interaction is independent of 
spin (i.e. when the scattering lengths are equal). It is similar to the molecular field of a magnetically ordered system 
and plays a central role in the collective spin dynamics. At T = 0, however, when the noncondensate is absent, this 
term vanishes. The effect of different scattering lengths is to modify the longitudinal component of the external field 



n'^ ^ (r?^, ny, -i- 2A„ + aj. 



(28) 



The mean field frequency shifts (or "clock" shifts) due to the condensate Ac and noncondensate A„ (which is multiplied 
by the Bose-enhanced factor of 2 in Eq. (|28|l Ql) are 



Ac = [5ll"-cl +5l2"c2 - (522nc2 + 91211 cl)]/h. 



(29) 



A„ = [giini + 512^2 - (g22n-2 + gi2hi)]/h, (30) 
Finally, the anomalous correlation appearing in Eq. 123|) of the noncondensate is defined by 

mjj(r,t) = {i;j{r,t)i;,{r,t)). (31) 

The three field operator average ('(/'j(r, t)-0j(r, i)V'i(r, t)) in Eq. includes higher order correlations that will give 
rise to a binary collision integral describing the exchange of atoms between the condensate and the noncondensate. 
Before evaluating it explicitly, we next consider the dynamics of the noncondensate. 

The equation of motion for the noncondensate field operator tpi can be derived directly from Eq. (|18|l and Eq. H23|l : 

j 

j 
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+ X! 9vi'4'ji'ji'i - 'njji'i - ntjipj), (32) 
j 

where the effective potential [/„(r,t) and coupling field r2„(r, i) (analogous to Uc and ftc for the condensate) include 
the mean-field interaction of the noncondensate with itself and with the condensate: 

Un = C/oxt + + 522"-2 + ^(»^1 + n2), (33) 



n^ = n'^ + ^ Tr£(np + n). (34) 

The modified external field including the mean-field frequency shifts (or clock shifts) due to different scattering lengths 
is given by 

= (^!",^]^^^" + 2A„ + 2A,), (35) 

where Ac and A„ were defined in Eq. (|29|) and Eq. H3U|) . We note that this modified field for the noncondensate 
differs from Eg. H28|l for the condensate by the factor of 2 in front of Ac. These clock shifts were measured in recent 
experiments j4l|. We also note that the second term in Eq. (|34|l differs from that found for the condensate. In Eq. I|25|l 
for the condensate, the noncondensate spin couples to the condensate due to the exchange term in the mean field 
interaction, but the condensate itself does not contribute to this effect. In contrast, in Eq. (|34l) both the condensate 
and noncondensate contribute to this term. Above Tc, when the condensate is absent, this spin mean field term is 
still present and plays a dominant role in the collective spin dynamics of the noncondensate (it is in fact responsible 
for transverse spin waves in the thermal cloud [13, 113). 

It is useful to introduce a time evolution operator 'U(t,to) (distinct from \L{t,to) given in Eq. 114() ) that evolves ipi 
in time according to Eq H32|l 

^P^{r, t) = V){t, to)M^, to)U{t, to). (36) 
This operator Hit, to) evolves according to the equation of motion 

in^Uit,to) = He«U{t,to), (37) 
at 

where the effective Hamiltonian Hcs is defined as 

Hcs = Hq + H' , (38) 



(39) 



In writing Eq H38(l . we have separated the Hamiltonian into a part that describes the mean-field dynamics ifo and 
the remaining part iJ', which can be treated perturbatively. The mean field part is given by 



-— v2 + c/„ i + -f7„-o: 

2m / - 2 — 



J' 



where the effective potential Un and coupling field Vln have been defined in Eqs. (|33|I - H35|1 . 
The four terms contributing to H' are 

H[ = -^gtj / dr i^^iUrj 4- ^Jn^J + ^*m^j -f (i/'jV'jV'i)) 'ip] + H.c. 

ij 



(40) 



(41) 
(42) 
(43) 
(44) 
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In the above formulas, the arguments of the field operators are (r, io), i-e., ip = ^o)- Note that the effective 
Hamiltonian HcS has time dependence through nij{r, t), $i(r, t), and mij(r, t) (apart from the possible time variation 
of the external potentials), which we have suppressed so as not to confuse these quantities with the Heisenberg 
operators. We emphasize here that we have made no additional approximations in introducing U, that is, one can 
show that Eqs. 138l) - (|44|l . together with Eqs. (|36|l and l|37() . reproduce the original Heisenberg equation of motion in 
Eq. 123. 

For an operator 0{t) constructed from the noncondensate field operators ijji{r, t) and 'i/'|(r, i), one can show that 

tr/5(to)0(i) = iTp{to)U\tM)0{to)U(tM) 
= iT~p{tM)0{to) 

^ {0)t. (45) 
This defines the effective noncondensate density operator 

p{t,to)=n{t,to)p{to)U\t,to), (46) 

which obeys the equation of motion 

i^j^'pit^ to) = [Hc«{t),pit, to)]. (47) 

Having defined the time evolution of the noncondensate field operator, we now turn to the derivation of the kinetic 
equation for the noncondensate. We introduce the Wigner operator 

W,,{r,p,t)= ydrVP-^'/ft^t(j. + r72,t)^,(r-r72,t), (48) 

and define the semiclassical distribution function 

l^„(r,p,t) = (z|^(r,p,i)|j) =trp(to)T^y(r,p,i) =tr/5(i,io)M^y(r,p,to). (49) 

Knowledge of this function allows one to calculate various noncquilibrium expectation values, such as the nonconden- 
sate density introduced in Eq. lf?7|) 



n^,{r,t)= I ^^W,,{r,p,t). (50) 



The equation of motion for Wij given from Eq. H47|l is 



^W,jir,p,t) = ^trp(t,io)[VFy(r,p,io),ffcff(i)] 

= ^trp{t,to)[W^j{r,p,to),Ho{t)] 

+ 4tr/5(t, to)[W,, (r, p, to),H'{t)]. (51) 
in 

The first term on the right hand side defines the free-streaming operator in the kinetic equation. With the assumption 
that U„{r,t) and ri„(r, i) vary slowly in space, then we can neglect the quantum corrections to the free-streaming 
term J9j , which leads to the semi-classical result 



Vrffi- 7;{V^£,Vp^}- =1, (52) 



dt m —2' '- h' 

where [ , ] and { , } on the left-hand side represent the commutator and anticommutator for the 2x2 matrices. Here 
we have defined the 2x2 matrix for the noncondensate effective potential 

2(r, t) = t/„(r, t)l + |d„(r, t) ■ £. (53) 
The term on the right-hand side is given by 

= ^trp{t,to)[mj{v,p,to),H'it)]. (54) 
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To summarize the derivation so far, we now have the condensate equation of motion Eq. (|23|l and the kinetic 
equation for the noncondensate Eq. (|52|l . Both of these equations depend on higher order correlations, which enter 
through the three field average (V'J(r, t)'0j(r, i)?/'i(r, f)) in Eq. H23|l and through the collision term / in Eq. (|52|l . In 



order to obtain a closed set of equations for ^i{r,t) and ^(r, p, i), we must truncate the description by making a 

Markovian ansatz about collisions. This can be carried out using the machinery of perturbation theory. 

To carry out perturbation theory, we first write down the formal solution to the equation of motion Eq (|47|) for 



dt%{t,t')[H'{t'),p{t',to)]Ul{t,t'), 



where 'Uo(i,to) is the unperturbed evolution operator 

Uo{t,tQ) = Texp 



dt'Hoit') 



to 



(55) 



(56) 



with T being the time-ordering operator. Iterating Eq. 155() to first order in H' , one has 

p{t,to) ^ Uo{t,^to)p{to)Ul{t,to) 

- [ dt%{t,t')[H\t')Mt\to)Kto)Ul{t' ,to)]Ul{t,t'), 

" J to 

Using this in Eq. 145|) , a statistical average of an operator O constructed from ipi can be expressed as 

{d)t = tvp{to)[ul{t,to)d{to)Uo{t,to) 



(57) 



to 



d^'uS(^^^o)[U^,(^,^')0(^o)1io(^,^'),-?^'(^')]1io(^^^o)}, 



(58) 



We may now use this perturbative solution to calculate the quantities {■4'^j{r,t)ijjj{r,t)ipi{r,t)) and I, which involve 

averages of products of three and four field operators. As in ZNG [s^, we effectively calculate these collision integrals 
to second order in g^ , while explicitly keeping interaction effects in excitation energy and chemical potentials only 
to first order in gij. Within this approximation, one can consistently neglect the anomalous correlation function, i.e., 
we set ifiij = 0, and thus there is no contribution from H2{t). A detailed discussion on the anomalous correlation 
function in a single-component BEG is found in ZNG j^^] . 

In Appendix^ we provide a detailed derivation showing how the three and four field operator averages are reduced 
to binary collision integrals using the perturbative solution in Eq. H58|) . There, we carry out this calculation for the 
general case of unequal scattering lengths. However, for simplicity, and motivated by the near degeneracy of the 
scattering lengths in ®^Rb, here we focus on the case an = 022 = 012 = a (with g = AT:Ti^a/m). In Section llVI we 
consider how having different scattering lengths modifies the kinetic equations. 

The dynamics of the condensate and noncondensate are described by equations of motion that are coupled through 
mean field interaction terms and collision integrals that describe the exchange of atoms between the condensate 
and noncondensate. The condensate spinor $(r, t) is described by a generalized finite temperature Gross-Pitaevskii 
equation 



2 m 



{Kc - iR)^ 



1+2"'= 



(59) 



and the noncondensate atoms are described by the kinetic equation 



dm I 



dt 



Here we have defined the condensate Hamiltonian _H_c and the semiclassical noncondensate Hamiltonian H r, as 



Sc(r,t) = 



2m 



Uc{r,t) 



(61) 
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l+-f1„(r,t)-£. (62) 

The dissipative term in the GP equation Eq. (|59|l and the two coUision integrals in the kinetic equation Eq. JHOJ, 
for the case of equal scattering lengths, can be obtained from the more general expressions Eqs. (|B47IB49p in Appendix 

m 

xS{pc + Pl - P2 - Ps) 

X {(1 + ffi3)Tr[m(l + W2)] - ^Tr[(l + W9 ] 

+ (1+ ^2)^(1 + ^3) -^(1+^)^3}, (63) 



Rnir,p,t) = 



2m 



+ Un{r,t) 



xS{p + P2 - P3 - P4) 

X ({1 + m, m}Tr[(l + ^2)^3] 
- 1 + m }Tr [^2 ( 1 + ^3 ) ] 
+{1 + m Ksil + K2)K4}- 

{ffi,(l + ^3)^2(1 + (64) 



LlWji = ^ J J^niif J '^^'^ J dp36{ec + ip, - ip^ - ip^) 

xS{pc + Pl - P2 - P3){['^(P - Pl) - S{p - P2)] 

X [{1 + ^2}Tr(^3^) - {m, 1 + E2}Tr[(l + ^3)^] 
+ [(1 + )W-,n,. Wp + E2ncffi3(l + mi) 
-mil + K3)Rcil + K2) - (1 + K2)^il + K3)Kl] 
-S{p - P3) (fe, E3}Tr[(l + m)ffi2] 

-{Up, 1 + E3}Tr[m(l + K2)] 

+nrW9 (I + W, W3 + W,3 (1 + m )W_',nr 

-np(l + E2)m(l+E3) - (1 + ^3)^1(1 + ^2)1^)}, (65) 



where ^ = W,v, p, t) and ^ = ^(r, p^,t). 

The local condensate and noncondensate energies are defined by 



ec(r, t) = Re[$^^$]/nc, (66) 



ep(r,t) = Tr(^„£)/n, (67) 

while the condensate momentum is defined by 

p, EE TOV, EE V[1>* V<I>. - $,V<I>*]. (68) 

Ziric ^ — ' 

i 

More explicit expressions of the above three quantities are given in Eqs. H73|) . ljB17|l and H72|) below. 

At this stage we pause to discuss various cases of the spin- 1/2 kinetic equations H59|) and H6U|) and make contact 
with some earlier work. We first remark that in the case of T > Tc, the above result for X„ agrees with the collision 

integral derived by Jeon and MuUin [l^f, and thus our kinetic theory above Tc reduces to the Jeon-MuUin theory. 
Furthermore, if one considers temperatures well above Tc where the Bose enhancement of collisions can be neglected, 
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one recovers the kinetic equation for a classical spin-1/2 gas introduced by Lhuillier and Laloe (we do not include 
the spin rotation term in the collision integral, as we discuss in Appendix 

If we set the internal coherence of the noncondensate W12 = to zero, and take the coupling field to zero (l — 0, 
then our theory becomes equivalent to that for a binary mixture of two different atomic species (with the same mass). 
For the case where the population of one of the components is zero, our kinetic theory reduces to the ZNG kinetic 
theory for the single-component Bose-condensed gas [231 • 

In this paper we are interested in the general case where the internal coherence of the noncondensate is finite. 
Comparing Eqs. (|59|l and (|60|l with (3) and (4) of Jackson and Zaremba we see that the structure of the spin-1/2 
kinetic theory is the same as the ZNG theory for a single component gas. In particular, the collision integrals Lc 

and /„ are generalizations of the C12 and C22 collision integrals given by Eqs. (23a) and (23b) in ZNG "35^. The 

spin-1/2 Bose-condensed gas, however, is a much richer system due to the extra spin degrees of freedom. For a given 
excitation symmetry, the number of collective modes for the spin-1/2 gas is three times as large compared to the 
single component gas, for one has total density, relative density (or longitudinal spin), and transverse spin dynamics 
to consider. Similarly to the single component gas, one must also consider the relative motion of the condensate and 
noncondensate for each of these degrees of freedom. In order to better understand the various types of motion, it is 
extremely useful to recast the theory into a form that separates out the spin and total density degrees explicitly. 



III. SPIN HYDRODYNAMICS 



The condensate spinor order parameter can be described in terms of density and spin variables, which obey equations 
of motion exactly equivalent to the spinor GP equation. The condensate spinor can be written as |39| 

- [^lirU) ) = ^J^y^^'-'' [ e'^=(-*)/^sin%i) ) ' ^''^ 

where ric = ?ici + «.c2 = I'&iP + 1*5*21^ is the total condensate density, and ac is the overall phase of the condensate 
wavefunction. The two phase angles (pc and 9c define the spin (Bloch vector) of the condensate 

5e(r,i) = Tr{£^e(r,i)} = n,(r, i)Me(r, i), (70) 

where is the reduced condensate spin vector 

Mc — (cos0csin6'c, sin (/)c sin 6'c , cos^c)- (71) 
The local velocity flow Vc(r, i) for the total condensate density nc(r,t) is given by (see also Eq. (|68fl ) 

Vc = — [Vac - ]- cos6'cV0c ) . (72) 



We remark that the velocity field of a spin-1/2 order parameter is not irrotational (i.e. V x Vc 7^ 0) in general, a well 
known result. Finally, the local energy of the condensate is given by (see also Eq. (IGGfl ) 

e,(r,t)^^,(r,i) + ^^^^^, (73) 
where the explicit formula of the local chemical potential for the condensate is given from Eq. (PTT|) : 

^ _ hy^n {r t) ^ ^ ^ ^ ^[VA4(r,t)]^ (74) 

The equations of motion for the condensate total density, velocity, and spin density can be obtained directly from 
the spinor GP equation IjSQI) using Eq. 1)691) 



— + V.(n.v.)= — 



(75) 

coll 
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dt 



dec 
dx,, 



— Q . t 

2 " dx„ 



n dMr 



dMr d 



dMr, 



Amric dx^ dxi, \ dxi, 



dMr 
\^ dx^, 

1 / - dMr 

5 "-"is: 



Tr(££), 



(76) 



at m h ot 



(77) 



coll 



where and u in Eq. (76) are Cartesian components and repeated subscripts are summed. Here the noncondensate 
spin is defined as 



(78) 



The condensate spin current is given exphcitly in terms of the other condensate variables by 



Jc = mScVc - - — Sc X VSc 
2nc 



(79) 



We note that Jc(r, t) is a tensor field, and the operation V • Jc(r, t) should be interpreted as 

d d d 

V • Jc(r,i) = —Jr^{r,t) + —Jry{Y,t) + —Jrz{r,t). 

The collisional contributions in Eqs. (|75|) and H77|) are given by 



dur 
'dt 



coll 



= --Tr{nc,^}, 



(80) 



(81) 



dSr 

'dt 



-Trg{^,2}] 



(82) 



We now turn our attention to the noncondensate kinetic equation. The noncondensate hydrodynamic equations 
can be obtained by taking moments of the kinetic equation, and inserting a local equilibrium form for ^(r, p, t). The 

first three moment equations define the spin density, spin current, and spin pressure of the noncondensate component: 



Sn{r,t) 



dp 



:TTaW{r,p,t), 



(83) 



J«(r,t) 



dp 



rpTr£K(r,p,i), 



(84) 



P^..{r,t) = 



dp Pf^iPu 



Trff ^(r,p,i). 



(27r?i)3 TO 

It is straightforward to derive the general moment equations from the kinetic equation Eq. H60|) : 

at m Ti at 

coll 



(85) 



(86) 



9J„„, dP„ 



dU„ 



dt dxi, dx. 



h _ di}„ dJ„ 



2 dxfj, dt 



coll 



(87) 
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The subscript /i in J„p signifies a coordinate space vector component of J„. Here tlie collisional contributions are 
given by 



dt 



coll 



d3n 

dt 



coll 



(89) 



We note that the equation for Sn in Eq. (|86|l involves the colhsional contribution described by the following collision 
matrices 



r„ 



dp 



{2'Khy 



Mm 



(90) 



dp 



Mm- 



(91) 



Eq. H88|l can be written in terms of these collisional matrices as dSn/ dt\co\\ ~ TrCT(£c the case of equal 

scattering lengths we are considering here, one immediately finds ^„ — 0, i.e., collisions between noncondcnsate 

atoms conserve the density and spin density. In contrast, collisions between the condensate and noncondensate do 
not conserve the density and spin density of each component since atoms can scatter into and out of the condensate. 
One can show that the dissipative term appearing in the GP equation is related to Fc, and thus to /c, according 

to 



1 



(92) 



coll 



where dnc/dt\co\\ is the collisional contribution to dn^/dt. The relationship between /c and Fc in Eq. (|91|l . along 

with Eq. (|92|l . ensures the conservation of total spin and atom number. The collision matrix Tj., which is a source of 

damping of density and spin excitations in the condensate, vanishes only when the condensate and noncondensate are 
in complete local equilibrium with each other, as we will show later. When we consider unequal scattering lengths e.g., 
5ii 7^ 5i2, En also makes a finite contribution to Eq. (|86|l . giving rise to relaxation of the transverse spin component. 

We discuss this in Section Hvl 

In order to reduce the general moment equations to a closed set of equations, we assume that rapid X„ collisions 

brings the noncondensate into local equilibrium. To be in this collision-dominated hydrodynamic regime, the collective 
mode of frequency uj must satisfy lut < 1, where r is some appropriate relaxation time for reaching local equilibrium. 
One can show that the appropriate local equilibrium distribution function is given by |18| | 



^"=^(r, p, f) = ^ i [1 + se„(r, t) ■ ^/.(r, p, t), 



(93) 



where s is +1 or —1 and e„ is a unit vector representing the local direction of the noncondensate spin moment. One 
can express e„ using the transverse and longitudinal angles of the spin vector as 

e„ — (cos (j>n sin , sin 0„ sin 0„ , cos 0„ ) . (94) 

The distributions for the dressed states (defined along the local polarization axis e„) /±(r, p, t) are given by 



2m 2 



1 



(95) 



Here jigir^t) is a chemical potential for a dressed state s, which may vary in position and time. With this local 
equilibrium distribution, the local spin density is given by 



Sn = e„(n+ - n_) = e„nM„, M„ = (n+ - n_)/(n+ + ri_). 



(96) 
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where 



Arp 



exp 



(97) 



and Xt = {^-Kh? /mk^TY^'^ is the thermal de Broghe wavelength. For later use, we also express the condensate 
component in the analogous form 



^(r,t) 



[l+ec(r,t) -ol. 



(98) 



We note that the condensate spin is always fully polarized with the local direction Cc, and thus one always has Cc = Mc- 
It is straightforward to verify that this local equilibrium distribution satisfies L„ [M^^"*^] = 0, independent of fig and 

e„ (this condition indeed defines the local equilibrium solution). In contrast, Xc in general does not vanish, i.e., 



dp2 / dp 



2n J (27r?i)3 

Xi5(mVc + Pl - P2 - P3)5(£c + Epi - - Spa) 
X _ g-/3[A,/-™(v„-Ve)V2-Pc]| 

s,s' = ±l 

x{[J(p- Pl) - ^(p - P2)](l+ se„ •£)(! + s'e„ • e;)(l H 
x[(5(p - Pl) - 6{p - P3)](l - Ss,s')a, - [ec - e„(e„ • Ec)] 
-S{p - P3)[(l + s'e„ • ec)l+ (cc + s'e„) -£1(1 + Ss^s')} 

X(l + /ls)/2s/3s' 



(99) 



This equation is analogous to Eq. (40) of ZNG [33 for the C12 collision integral in a single-component Bose gas. If 
we impose the condition that the first square bracket vanish for both components s = ±1, one must have = 
+ Tnixn — ^cY I'i- This condition is too restrictive, since it results in f+ = f-, which occurs when the polarization 
of the noncondensate vanishes. In order to satisfy Lc — while allowing for a finite noncondensate spin moment, the 

condensate and noncondensate local spin polarizations must be parallel e„ = Cc (or antiparallel e„ — — ec). In this 
case, one can show that all the contributions from the /_ (or /+ in the antiparallel case) component vanish. Then 
local equilibrium between the condensate and noncondensate is specified by the conditions (for the parallel case) 



(100) 



/X+ = + y (Vn - Vc) . (101) 

The role of the condition for the chemical potential in Eq. I|101() has been discussed extensively for the single- 
component Bose gas [ssl [tJ ; for example, coUisional equilibration between the condensate and noncondensate gives 
rise to damping of condensate collective modes [s^ • 

In the spin-1/2 case, there is an additional condition in Eq. (|100|l for the the spin orientations. This additional 
condition means that there is a coUisional relaxation process related to the relative orientation of the condensate and 
noncondensate spins, which will give rise to damping of condensate spin dynamics. This relaxation process comes 
into the dynamical equation of motion for the condensate and noncondensate spins [see Eqs. (|88|l and (|92|) ] through 

X(5(mVc + Pl - P2 - P3)^(ec + £pi - £p2 ^ ^Pa) 
X ^ _ g-/3[/i,,-m(v„-v,)V2-M<:]| 

s,s'=±l 

[(1 + s'cn ■ ec)l + (ec + sen) -£1(1 + Ss,s') 

X(l + /l.)/2s/3s'. (102) 

The role of this term in the spin dynamics will be more explicitly seen later. 

To obtain the equation for the spin current, we must include the deviation from local equilibrium. For simplicity, 
we consider the situation where the condensate and noncondensate components have no particle currents. We then 
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follow the linearization procedure given by Jeon and Mullin [l^ . expanding the distribution function around local 
equilibrium as 

^ = ^I'^q + SK- (103) 
The lowest order correction is found to take the following variational form p^ : 

where ^|:'^ is perpendicular to e„, i.e., ^/J"-^ • e„ ~ 0. The functions T/ijl and ip'^ are associated with the spin current 
through 

where ji and denote the longitudinal and transverse components 

J„ = 4 + Jn = en Jll + Jn , 4 ^ e„ , (106) 

and Ps is defined by 



rfp r ksT 

With this approximate form of the distribution function, the spin pressure tensor becomes diagonal: 

P^., = en6^,Y.''^- = ^^^^^^^ P = &„EsP.. (108) 

s s 

We note that P and Sn are related through Eqs. (|96|l . H97(l . H107(l . and p()8(l . along with the condition hi{r,t) + 
n2(r, t) = no{r). While Jeon and Mullin |13i found this variational solution for a noncondensed gas above Tc, it is also 
applicable to the thermal gas in the Bose-condensed phase, because our kinetic equation in Eq. H6Q|I is formally the 
same as that of Ref. (i.e., the condensate contributions are buried in the effective potential ^^ and the collision 

integral). We then obtain the following equation for the noncondensate spin current: 



-f + VP + vuA + -vn^n - n„ X J„ = -f 

ot 2 ot 



(109) 

coU 



To evaluate the collisional contributions in the hydrodynamic equations, we must linearize the collision integrals in 

um = ici^'"'^] + 5Lc[5m^ (110) 



■^l and "0 



(111) 



The contribution from ^'""^ in Eq. piOfl leads to L given in Eq. (|l()2|l . In [O we give detailed calculations of 

the linearized collision integrals. Here we simply present the main results. Assuming that the condensate spin and 
noncondensate spin are almost in local equilibrium with each other, we linearize Eq. p02|l in ^diff = ~ A*c and 
e„ — Cf. to obtain 



dSn 

dt 



coU 
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We also obtain 



dJn 

dt 



coU 



dp 



rpTr 



£((5/„ + 6Lc) 



n 



•D 



(113) 
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Explicit formulas of the four relaxation times f| , , rf^, and are given in Appendix [CI Here we note that the 

spin diffusion relaxation times rf^ and involve the contribution from Xc term, which represents collisions between 

the condensate and noncondensate. 

In summary, we have obtained the following coupled spin hydrodynamic equations, assuming no overall mass 
currents for the condensate and noncondensate: 



1 



n' X Sr 



9l2 



Sn X Sr 



(114) 



2nr 



(115) 



dt 



-V • J„ 

m 



9l2 



n 



Cc), 



(116) 



at 



VP + VUnSn + - VSl„n - fl„ X 



D 



'D 



The modified external field can be written as 

1 



^'c = ^ + ^^[(ffii - ff22)(ne + 2n) + (gii +522 - 25i2)(^," + 25^)], 



(117) 



(118) 



n'^ = n + z-[{gn - g22){nc + n) + {gn + 922 - 2ffi2)(5e + S^)], (119) 
while the effective field for the noncondensate is 

(in = ^'n + ^{Sc + S^). (120) 

Ti 

Here we remark on some important features of the spin hydrodynamic equations. 

1. The condensate spin current involves the spatial gradient of the condensate spin density itself. This term exert 
a torque on the condensate spins, which causes spin untwisting observed in the JILA experiment [39j |. It is 
important to note that this term dose not involve the interaction, and is a sort of quantum pressure term (also 
referred to as the "quantum torque" HJl). This is in sharp contrast with the noncondensate spin dynamics. The 
exchange mean field plays a crucial role in the noncondensate spin current, which causes collective spin waves. 

2. The above equations involve coupling between the condensate and noncondensate spins through the exchange 
mean- field terms, such as {gi2/K)Sc x Sn, which exerts a mutual torque when the condensate and noncondensate 
spins are not aligned. 

3. The two components are also coupled through the coUisional relaxation terms. These terms try to bring the two 
components into complete local equilibrium i.e., /i+ = /^c and e„ = Sc- These conditions determine the length 
and the direction of the spins. 
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IV. EFFECT OF DIFFERENT SCATTERING LENGTHS 



The collision integrals given in Eqs. (|63I65() were obtained for the case of equal scattering lengths gij = g, corre- 
sponding to SU(2)-symmetric interactions. For this special case, collisions make no explicit contribution to the decay 
of the net transverse spin. In this section we consider the effect of unequal scattering lengths, which turn out to 
contribute to the decay of the net transverse spin. 

We first consider the role of /„ integral. From Eq. (|90|l . the collision matrix r„ is given by 



TT f dpi f dp2 f dp3 



dp4 



h J {2TTh)^ J {2Trh)^ J {2Trhy 
xS{ep, + Ep., - ip^ - SpJS{pi + P2 - P3 - P4) 

x[{l + ffipi),A<(P2;P3,P4)} 

-{ffipi),4^(p2;P3,P4)}], (121) 

where the explicit formulas for the matrices A-^ are given by Eq. (jB31|) . After some rearrangement, we obtain 

_ n f dpi f dp2 f dp3 f 
- hj {2nhrJ {2nnfJ {2^%^ J '"'^^ 

x5{sp^ + ip^ - £p^ - ep4)i5(pi + P2 - P3 - P4) 



klm 

-W,< (pi)iy< {p2)W>ip,)W>^ (P4) 
+ W,> (Pl ) (p2)T^<fc (P3) W^,< (P4) 



-W<MW^iiP2)W>MW>{p,)] } (122) 

One immediately finds that Tnu = regardless of the form of ^(p). This means that the collisions between 

noncondensate atoms conserve the number of noncondensate atoms in each hyperfine state. We now focus on the 
finite contribution of the transverse component, namely 

TT f dpi f dp2 f dps f 



n J {2TThf J {2TThf J {2TThf 

xS{ep, + ip^ - ip^ - epJS{pi + P2 - Ps - P4) 

:{E(5fc2 - 9lk)gml[WiUPl)W^liP2)W<{p3)W<2{p4) 



klm 

-W<„{Pl)W,<i{p2)W>{p,)W>2{pi) 

+w>^ {pi)w> {P2)w<MWi'i (P4) 

-W<,{pi)W<{p2)WlMW^2iPm (123) 

The above formula vanishes if gn — gi2 — .922, i-e., the transverse spin is conserved during collisions. If the scattering 
length are not equal, collisions do not conserve the transverse spin. 

Assuming a small transverse spin component, i.e. \Wi2\ <C Wn, W22, we expand r„i2 in W12 to first order 



27r f dpi f dp2 f dp. 



dp4 



h J (27r?i)3 J (27r?i)3 J {2T:Kf 

X6{ep, + £p^ - - £pi)5{pi + P2 - P3 - P4)W^12(Pl) 

X 5l2(.9l2 - gil){[l + /1(P2)]/1(P3)/2(P4) 

-/l(P2)[l + /l(P3)][l + /2(P4)]} 
+922(922 - gi2){[l + /2(P2)]/2(P3)/2(P4) 
-/2(P2)[1 + /2(P3)][1 + /2(P4)]} 
+912(912 - .922){[1 + /2(P2)]/2(P3)/l(P4) 
-/2(P2)[l + /2(P3)][l + /l(P4)]} 
+511 (ffll - 5l2){ [1 + fl (P2)]/l (P3)/l (P4) 

-/l(P2)[l + /l(P3)][l + /l(P4)]}l, (124) 
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where fi = Wu. 

To obtain a closed form in terms of the transverse spin component, we assume the local equilibrium distribution 
given in Eqs. 1)93 11 - 1)95(1 . Then /i, /2 and W12 are written in terms of /+ and /_ as 



/i = ^ [(/+ + /-)+ cos 0„ (/+-/_)], 



(125) 



/2=2[(/+ + /-)-cos0„(/+ -/-)], 



(126) 



W12 = -(/+-/-) sin 0„e 



With this form of W12 , we find 



dt 



coll 



ni2 

T2 



(127) 



(128) 



In order to obtain the explicit form of the relaxation time T2, we consider the following two cases. First, we assume 
that the spin direction is almost along the z direction i.e., 9n ~ 0. To first order in 0„, one has 



/2«/-, W,2 ^ \{fl - f2)ene-'^" 



(129) 



The relaxation time T2, which is a kind of "T2" lifetime, is given by 



1 



27r 



dp I 



T2 [ni - h2)h J {2Trh)^ J (2nh)^ J (2nh) 



dp2 



dp3 



dp4 



^P3 - ^pJHpi + P2 - P3 - P4) 



(511 - 512)'{[1 + /2(P1)][1 + /1(P2)]/1(P3)/1(P4) 

-/2(Pl)/l(P2)[l + /l(P3)][l + ,/l(P4)]} 

+ (522 - 512)'{[1 + /2(P1)][1 + /2(P2)]/2(P3)/1(P4) 

-/2(Pl)/2(P2)[l + /2(P3)][1 + /l(P4)]}] . (130) 

If the spin is almost polarized along the longitudinal direction, i.e. /2 ~ 0, the above relaxation time reduces to 



1 



27r 



dp I 



T2 nih J (27r?i)3 J (27r?i)3 J {2i:K) 



dp2 



dp3 



dpi 



£p4)'^(Pl +P2 - P3 - P4) 



x(gii - .gi2)^[l + /i(P2)]/i(P3)/i(P4). 
In the opposite limit of the small longitudinal polarization /i w /2 and d„ ~ 7r/2, we write 

TT 

0n = -+59n, fi+=il + Sfl, fl^=fl-6fi. 

To first order in the small variables, one has 



with 



We then obtain 



1 



1 

T2 



93/2(2) An f dpi f dp2 f dp3 



(131) 



(132) 



(133) 



(134) 



dpA 



9i/2{z)nhJ {2TThY J {2nh)3 J {2-KhY 
xS{£p, + ip., - £p^ - ep^)5{pi + P2 - P3 - P4) 
X [(511 - 512)^ + (522 - 512)^] 
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X[l + /(pi)][l + /(p2)]/(p3)/(p4), 



(135) 



where h — ni + n2 ^ (2/ Af-) (73/2(2). One can show that Eq. H13()(l smoothly goes over to Eq. H135(l . if wc take the 
hmit f« /2_ ~ fi. 

We now consider the simplest case of the high-temperature limit T 3> Tc- In this limit, the local equilibrium 
distribution takes the Maxwell-Boltzmann form 



m'^^ir^p^t) « i[l + M„(r,t) •a/MB(r,p,i), 



with 



/MB(r,P,i) = exp ~I3 



rp^ 

2m 



/i(r,i) 



(136) 



(137) 



Assuming /mb ^ 1 and keeping the quadratic terms in Eq. H123|) , we obtain the following simple form of the relaxation 
time 



1 

T2 



dpi 



dp2 



dp3 



rdp4 



nhj^ {2TThf J {2TThf J {2TThf 

X(5(£pi + - - £p4)'5(Pi + P2 - Ps - P4) 
x{[(5ii - 5i2)'(l + M„-) + (522 - 5i2)'(l - M„-)] 

x/mb(P3)/mb(P4)}, 



(138) 



where n = (1/A|,)e^^. This expression for the relaxation time is formally identical to that of the classical collision 
time, and thus one explicitly finds 



T2 



ncTcSVth, 



where Wth = \/ ?>k-QT / nm is the thermal velocity and (Jcff is the effective collisional cross section 

aeff = 2^[(aii - ai2)2(l + i\/,^) + (022 - ai2)'(l - M^)], 



(139) 



(140) 



rather than the usual ctcoU = Stto . One finds that this relaxation time is much longer than the mean collision time 



T2/tc - (a/Aa)^ 



(141) 



where Aa is the difference in the scattering lengths. For ^^Rb, one finds T2/tc ~ 1000. In the recent JILA experiment 
[2^. this gives T2 ^ 10s. 

We note that if \gii — 912! — I522 — 512 1, which is the case of ®^Rb, T2 in the classical limit is independent of 
the longitudinal polarization M^. However, in the quantum degenerate regime, we see from Eq. H130|l that T2 may 
strongly depend on M^. This effect might be observable in experiments near Tc- 

We next briefly discuss the case of below Tc- Since the /„ collision integral has no explicit dependence on the 



condensate, E„ is the same as discussed above for T > Tc. For Tc, integrating Eq. (jB39|l over momentum, we find 

rc = r(i'+r(2\ (142) 



where 



and 



p(2) ^ 
=cii 



dpi 



n 

xS{e 

X 
mkl 



dp2 



2TTh)^ J {2TTh) 



dp3 



(143) 



coll 



-Pl 



^P2 - £p3)Sim-Vc + Pl - P2 - P3) 



^9mi{gkj - gik) 
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+w>Mw<Mw<{p3) - w,< (pi)w^>,(p2)w->(P3)] 

+ncik[W>^iPi)W<{p2)W<,iP3) - W<^ipi)W>{p2)W<,{p3) 

+W>^{Pi)W<{p2)W<M - W^<,(Pi)W,>(p2)M^<„(p3)] }. (144) 

(2) 

One immediately finds that T],^^ = 0. This means that collisions between the condensate and noncondensate conserve 
the total density n = ni + n2 and the longitudinal spin = ni — 712, but do not conserve the transverse component 
of the total spin. Further explicit calculations of L in the unequal scattering length case is quite complicated, and 

thus we do not make such analysis in the present paper. 



V. CONCLUSIONS 



In this paper, we derived a set of equations which describes the dynamics of a trapped Bose-condensed gas with spin- 
1/2 internal degrees of freedom at finite temperatures. These equations consist of a generalized Gross-Pitaevskii (GP) 
equation for the spinor condensate order parameter $(r, t) , as given by Eq. H59I) , and a semiclassical kinetic equation 
for the noncondensate distribution function ffi.(r,p,t) in a 2 x 2 matrix form, as given by Eq. (|60|l . These coupled 

equations are the spin-1/2 generalization of the kinetic equations obtained by ZNG for a single component gas |35j| . 
which has been shown to accurately describe several different experiments at finite temperatures 72] . Both of these 
theories treat the thermal excitations semiclassically in the Hartree-Fock approximation. Our kinetic theory is also 
the natural extension to the Bose-condensed regime of the earlier theories developed by Jeon and MuUin 0, 0, |23| 
and Ruckenstein and Levy jl9l | to describe spin waves in dilute quantum degenerate gases. 

As a further development of our theory, from Eq. (|59|l and Eq. 160|l we derived a closed set of spin hydrodynamic 
equations for the two spin components in the absence of any overall mass currents, i.e. v,, = v„ = 0, as summarized 
in Eqs. I|114l) - (|117|l . These equations are based on the assumption that collisions between noncondensate atoms 
(described by the £„ collision integral) are sufficiently rapid to produce local equilibrium. Deviation from local 

equilibrium gives rise to the usual spin-diffusion relaxation in the noncondensate spin current, which now involves an 
additional contribution from the collisions of noncondensate atoms with atoms in the condensate (described by the Xc 

collision integral). Interestingly, the Xc collision integral also gives rise the mutual relaxation between the condensate 

and noncondensate spins, which suggests that the condensate and noncondensate spins have a tendency to become 
aligned. This is a new prediction of our theory. In addition to the collisional coupling between the spins of the 
condensate and noncondensate, the exchange mean field gives rise to a mutual torque coupling the two spins. This 
mean field effect was first pointed out by Oktel and Levitov 'tt'I. 

There are a wealth of different applications of the present kinetic theory to finite temperature dynamics of a spin- 
1/2 Bose-condensed gas. Perhaps the most obvious application is to study the collective spin- wave dyna mics below 
Tc, which was never treated in the earlier literature on spin waves in spin polarized gases [T^ ITM Il9l l2Clj| . From 
earlier studies of spin waves above T^, we know that the exchange mean field has a strong effect on the collective 
dynamics of the transverse spin in the thermal gas due to the molecular field contribution to f2„ — Q'^ + gi2Sn/'h. 
Below Tc, there is an additional contribution to the molecular field experienced by the noncondensate due to the 
condensate r2„ = -I- gi2{Sn + Sc)/h- The molecular field experienced by the condensate has a reciprocal coupling 
to the noncondensate fl^ = + 5i2'S',i/?i, but notice that the condensate by itself does not experience the exchange 
mean field: interactions do not lead to collective spin waves in the condensate. Instead, spin waves in the condensate 
will result from the kinetic energy cost to twist up the condensate, an effect that is encapsulated by the so-called 
"quantum torque" form of the condensate spin current = —{fi/2nc)Sc x VS'c- Together with the spin relaxation 
effects due to collisions, these various effects should lead to interesting new spin- wave dynamics below Tc 82]. 

Another important application, motivated by recent experimental observations j83| ] is to study the effect that the 
dephasing of the spins in the thermal cloud has on the spinor condensate after a 7r/2 pulse is applied. In the presence 
of spatial inhomogcneities in the external field Vl'^ir), the spin of the noncondensate can decay to zero S'n(r) — > due 
to the dephasing of individual spins of atoms as they travel with different trajectories in the trap. In contrast, all the 
atoms in the condensate are in the same spinor wavefunction and so it is always spin polarized in our theory. As the 
spins of the thermal gas dephase, the gas evolves from a spin coherent sample to a 50 — 50 spin mixture, which has 
a condensation temperature approximately 20% lower than that of the spin polarized gas ^92] . At fixed temperature, 
we then expect the condensate fraction to decrease as the noncondensate spin decoheres. This "condensate melting" 
|83| is a finite temperature effect that requires a kinetic theory description that treats the scattering of atoms between 
the condensate and noncondensate, and so is an ideal application of the present theory. A related question arising 
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from this study concerns the vahdity of the present kinetic theory for describing the transverse spin dynamics of the 
condensate when it is subjected to such strong phase fluctuations. Indeed, two recent studies of the ground state of 
a spin-1/2 Bose gas suggest that in certain hmits, the condensate may become fragmented which is beyond 

our treatment of the condensate in the present kinetic theory. 

A final proposed apphcation of the theory is to treat the dynamics of topological states, or spin textures, of the 
condensate at finite temperatures. Several recent theoretical studies have explored spin textures at zero tempera- 
ture [s^ Is^ . We note that the two-component vortex state proposed ^| and observed [s^ at JILA is also a spin 
texture. At finite temperatures, based on our kinetic theory, one expects interesting effects due to the mean field and 
coUisional coupling between the spins of the condensate and thermal gas. In particular, the same kind of mutual fric- 
tion t ype effects arising in a single component Bose-condensed gas when the condensate and noncondensate velocities 
differ [73 will also occur in a spin-1/2 gas. The generalization of this effect on spin dynamics should follow from our 
theory. With the kinetic theory we have presented, one should be able to explore the role of the thermal cloud in the 
creation and stability of topological excitations, as was done in the study of vortex formation in a single component 
gas Iz!]. 
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APPENDIX A: TWO-LEVEL SYSTEMS IN DILUTE ATOMIC GASES 

With a bosonic atom, the two- levels making up the spin-1/2 system will always be a particular subspace of a larger 
space of internal states. In the alkali atoms, this larger space is the set of hyperfine states, illustrated in Fig. 1. In 
order to be a viable candidate, the system should meet a few basic requirements: the two-levels should have very low 
loss rates and long coherence times; it should be straight-forward to couple the states with a radiation field; they must 
also experience the same confinement while still allowing for spin-dependent forces that can excite collective modes 
of the spin. 

There will always be collisions that take the system out of the two-level subspace, which we will consider as loss 
processes. Typically, the most severe source of loss comes from spin-exchange collisions 0, (dipolar and 

three-body losses, which should be examined as well, are typically much smaller). The best way to minimize this 
loss process is to choose two states for which the selection rules or energetics do not allow the colliding atoms to 
change their internal states to ones outside the subspace. In these spin changing collisions, the total molecular spin 
-Ftot and spin projection Mtot of the colliding pair are conserved. For example, for colliding atoms in the states (1, 1) 
and (1, —1), the total spin can add up to the values Ftot =0,1,2 and the spin projection is Mtot = 0. Although the 
selection rules allow for scattering into states in the upper F — 2 hyperfine manifold, these collisions are disallowed 
energetically at ultralow temperatures. However, the process (1, 1) -I- (1, —1) (1, 0) -I- (1, 0) is allowed, making these 
two states a poor spin-1/2 candidate in general. Based on similar arguments, one can show that spin exchange losses 
can be eliminated by using either of the two pairs of stretched states: (1, —1) and (2, —2) or (1, 1) and (2, 2), which 
can be coupled using a single photon microwave coupling. For ^^Rb, any pair of hyperfine states may be considered 
since the spin exchange loss rates are anomalously low. 

In order to confine the two states within the same magnetic trap, they must have the same magnetic moment. In 
this regard, the stretched states (1, —1) and (2, —2) or (1, 1) and (2, 2) would be poor choices in a magnetic trap. It 
seems that the states (1, —1) and (2, 1) of ^^Rb may be the only viable spin-1/2 system in a magnetic trap. In order 
to use the stretched states, an optical trap must be used, which exerts the same force on all of the hyperfine states. 
In general there are various possible schemes for generating spin-dependent forces; a detailed discussion is beyond the 
scope of this paper. 

Based on these considerations, we have narrowed the possible two-state systems to either pair of stretched states 
(1, —1) and (2, —2) or (1, 1) and (2^2) for the alkali atoms ^Li, ^^Na, *^Rb in an optical trap, or the states (1, —1) and 
(2, 1) of ^''Rb in a magnetic trap |43,|83. These are illustrated in Fig. 1. We note that at present, only the system in 
Fig. la has been explored in any detail in actual experiments, but the spin-1/2 systems in Fig. lb would allow for a 
wider range of parameters to be explored (such as the effects of largely different scattering lengths in •^^Na or having 
an attractive interaction in ^Li). 
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APPENDIX B: DERIVATION OF COLLISION INTEGRALS 

In this Appendix, we give detailed derivation of the collision integral £„ and Xc given in Eas. (|65|l and 164(1 . We 

closely follow the approach of KD [tJ and ZNG 35]. As usual, we assume that the macroscopic variables are slowly 
varying in space and time compared to the spatial and temporal scale of a collision event. We thus approximate the 
macroscopic variables near (r, t) as 

ny(r',i') « n^J{r,t), C/„(r',t') w [/„(r,t), 

nc(r',0«^c(r,t), Uc{r' ,t') ^ Uc{r,t) (Bl) 
For the condensate order parameter, we expand the three phases ac, ^c, to first order in space and time: 

a,{r', t') « a,(r, t) + -^{t' - t) + Va, ■ (r' - r) = a,(r, t) + 5a,, (B2) 



,(r', t') « </.,(r, t) + ^{t'-t) + V^, ■ (r' - r) ^ </.,(r, t) 



f)0 

,(r', t') « 0,(r, t) + -^(t'-t) + Ve, ■ (r' - r) = ^.(r, t) + 59, 



From Eq. Ht)9|l we then obtain 

$(r',t') « e'^ 

~ exp <y I 

— exp 

where 



1 — i^5(f>caf — i'^59c{cos4'cgy — sin(f)cg^) 



(Jacl - ^5<j)cSf - ^59c{cos(t)c^^ - sin(/)c£^) 
'^"4 - ^5<j)cMc ■ £cos 6'c - ^£ • (Mc X (5Mc) 



(B3) 
(B4) 



(B5) 



5M, = ^(i'-t) + VMe-(r'-r). 
ot 



One can also write Eq. HB5|I as 



(B6) 



^(r',<') « exp 



Pc-(r'-r)-£,(i'-i) 



^(r,t), 



where 



(B7) 



A/c cos 6*^ + Mc X VA/c 



(B8) 



^jdac 1 ^ 
I (9i = 2= 



M,^cos0, + Af, X — 



(B9) 



The spatial gradients of the phases arc related to the condensate velocity Vc through Eq. H68|l . Inserting the spinor 
form of the condensate wavcfunction into the generalized GP equation Eq. 159(1 . we also find 



t ,'dac 1 „ dcj>c 
M^--cos^.— 



-£c, 



(BIO) 



where the condensate energy is defined by (we neglect the anomalous correlation m^) 



ec(r, t) 
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2 8m 
+;r777T E 5yRek*('"' 

mv'^{r, t) 



(Bll) 



where iJLc{v,t) is the local condensate chemical potential. 
It is now convenient to introduce the Fourier transform 



(B12) 



where V is the volume of the system. For convenience on calculations, we work with the finite volume V and thus treat 
discrete momentum variable p. After deriving collision terms, we take the limit y — s- cx) to replace the summation over 
p with the momentum integral. With the local approximation in Eq. HB1|I . one can write the mean- field Hamiltonian 
Ho{t') as 



4 J P 



(B13) 



With this approximate form of Ho{t'), the free evolution operator becomes 



lXo(i, t') ~ exp 



(B14) 



In performing the perturbation calculation using Eq. I|58l) . one needs to evaluate '\i\{t, t')aipilo{t, t'). It can be formally 
written as 



Ul(t,t')a,pUQ(t,t') 



(B15) 



In principle, one can diagonalizc Eq. (|B13|) to write down an explicit solution for Eq. IjBlSp . However, doing so would 
be quite complicated. Following Jeon and MuUin (ISj . we use the simple approximation of introducing the local HF 
excitation energy ip, 



'U.l{t,t')aip'\lo{t,t') w exp 



-ep{t-t') 



(B16) 



where 



£p(r, t) = j- + Un{r, t) + -r!„(r, t) ■ M„(r, t). 



(B17) 



The above approximation corresponds to replacing s_ with Af„l. We also use the same level of approximation for the 
condensate order parameter: 



$(r',i') « exp <^ i 



^(i'-t) + ^Pc-(r'-r) 



(B15 



where the local condensate momentum is defined by Pc(r,<) = mvc(r,t). We can also express H'{t') in terms of the 
Fourier transform: 
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xgij [e'P--/'^e-'^=(*-*')/'^($*fi,j + ^*hji)aip + H.c, 

H2{t') ~ XI 51 ^Pi+P2,2pc 

1 ■P2 



ij P1.P2 



1/2 Opo+Pl 



P2+P3 



V Pi ,P2 ,P3 



^4(0 « 2l^E E 

ij Pl>P2,P3,P4 

-EE ^'J' i^ji"-]p°'iP + %J«lp«ip)- 



+P2 ,P3 +P4 Sij f^lpi ^jp2 ^jP3 ^iP4, 



(B19) 



(B20) 



(B21) 



(B22) 



y p 



As usual, we assume that there is no initial correlation so that we can use the Wick's theorem to factorize terms 
such as 



("lpi"]p2«fcp3"'P 



'4 /to 



(«lpi«fcp3)to(a]p2a;p4>to + (alpia(p4)to(a]p2a/cp3)to- 



(B23) 



As we noted before, the anomalous correlations are neglected. Within the present approximation of slow variation of 
the macroscopic quantities in space and time, we can use 



(«]pi«iP2)to « SpuP2Wij{r,p,t), 



(B24) 



in the calculation of the collision integral. 

We first consider the three-field correlation function that appears in the generalized GP equation: 

{i,l{r,t)ijj{r,t)Mr,t)) = X (alp,«,P2«fep3)te-^(P^-P=-P^>-/^ (B25) 

P1,P2,P3 

Clearly the first term in Eq. (58) makes no contribution to the three-field correlation function. We then find 

(°'lpi'^JP2 0fcp3)t 

d^'(^iS(^^ to)[uS(t, tOaipi«iP2afcP2iio(i, iO, 'ff'(^0]^to(^^ to))to 

_ __ / J^l/pi{£pi-Ep2-Sp3)(t-t')/n 

X (kp,«.P2«feP2,^;(i') + m')]ho- (B26) 

The time integral is performed by setting to — > — oo introducing the convergence factor e~''(*~* ), which yields 



: Trd{ec + ip-^ - ip2 - Sps) + iV 



Cq + Sp^ Sp2 Sp3 



(B27) 



Using the Wick's theorem, we finally obtain 

(«lpiajP2afeP3)t 



-K5{ec + ffpi - ip2 - ^3) + 



g-ipe-r/fi. 



£c + £pi £p2 ^P3 



PC+PI1P2+P3 



X gim^l{Wmi{Vl)[^jm + Wjm{V2Wkl + W«(P3)] 
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+ Wmi(.Pl)[5jl + Wjlip2Wkm + WkmiPs)] 
-[l5m» + Wmi{Pl)]Wjl{p2)Wkra{P3)} , 



(B28) 



where we omitted the arguments r and t for simphcity. We thus obtain the three-field correlation function 

I 



- y 

'y2 



Pl^P2,P3 



7r(5(£c 



'<5pc + Pl,P2+P3 y^9lm^l{W„,i{pi)[5jm + Wjm(p2)][4; + Wkl{p?i)] 



Im 



+ Wrm{Pl)[Sjl + Wjl{p2Wkrn + WkraiPs)] 
-[Srm + W^i{pi)]Wj,n{p2)Wkl{p3) 
-[S,ni + Wmi{Pl)]Wjl{p2)Wkra{P3)} , 

Using this result, one obtains the explicit formula for the three-field contribution in the GP equation as 

j 

■ I sr - ~ ~ ~ ^ -f 1 



(B29) 



Pl.P2,P3 

X'5pc+Pi,P2+P3 [A< (pi; P2, Ps) - A>(pi; P2, P3)]*i 



(B30) 



where we have defined the matrices 

4(Pi;P2'P3) = E-9fc^5.;{W^fc1(Pi)w|(P2)W,>(p3) 



kl 



with 



+ W<{pi)W>{p2)W>{ps)}, 



^<(p)=^(p), ^>(p) = l+^(p). 



(B31) 



(B32) 



This notation using > and < follows the Kadanoff-Baym formalism (actually, these matrices are directly related to 

the second-order self-energy in the Kadanoff-Baym formalism) l6ll Isof . We note that the matrix Af^ has the 
property 



A>(Pi;P2,P3) - [A>(pi;p3,P2)]* (B33) 

Eq. ljB30|l defines the matrix ^ that appears in the generalized GP equation Eq. I|59|l . 

The similar techniques and approximations are used to evaluate the collision terms in the kinetic equation. We first 
consider the contribution from H!^: 

-UTp{t,to)[W,,{r,p,t),H's{t)] 



P2+P3 



q Pl,P2,P3 



' {e'P"""/'' E[f''^-'^P.Pi-q/2*fcKp-q/2a«P2«fcp3)t 

k 

'9kjSp,p^+q/2^k{0'lpaip+q/2akp3)t - 3fei'5p,p3+q/2*i("Ip«fcP2a»p-|-q/2)t] 

-e*'''^''"^''E['^pi^p-q/25fcj*fc(4p3«lp2a«p+q/2>t 

k 

-^P2,P+q/2ff»fc*fc(a]p_q/24p3«»Pl)* 
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3+P4 



Finally, we consider the contribution from H'^: 

q \ fc pi,p2,P3,P4 

X [Spi ,p+q/2.9jfc («jp_q/2«Ip2 "fep3 "-ipi ) * 

^'^P4,p+q/2gkj{a]p^alp^akp3a^p+q/2)t]j 
- X!{5fej[^fejW^<fc(p) + ^fefeW^jj(p)] 

k 

-Qik [riikWkj (p) + fikkWij (p)] } . 
We now need to evaluate the four-field correlation function: 

("ipi«Ip2"fcP3««P4)t 

= tr/5(io) {t, ^o)a]pi alp^ akp^a.p^Uo {t, to) 
X ill {f, io) [a]p^ a^p^ afcp3 a,p, , {t')]Uo , to ) } • 



(B41) 



(B42) 

We find that the first term of order makes no contribution to Eq. IIB41|l . Using Wick's theorem to evaluate the 
second order term, we also find that mean field contribution in the last two terms of Eq. ljB41(l are canceled out with 
several terms in Eq. ljB42|l . The relevant terms are often referred to as "connected diagrams" . We thus obtain the 
relevant contribution: 



(«ipi'^lp2"fcP3"iP4)t 

i 

~ ~V 



-Pl ' ^P2 ^P3 ^P4 



-'Pl+P2,P3+P4 



X]5mn|VFmj(Pl)W„fe(p2)[(5fe« + VFfc„ (ps)] [(5,™ + Wi™(p4)] 

rnn 

+ Wmj{Pl)Wnk{p2)[Skm + WkmiPsWtn + Winipi)] 
-[Smj + W,nj{pi)][5nk + W^nfc (P2 )] W^fen (P3 ) (P4) 

-[Srnj + W,nj{pi)][5nk + W„fe (p2 )] W"fo„ (pa ) Wi„ (P4)}. 

Using the matrices defined above, we can write this H^^^ contribution to the kinetic equation as 

^trp{t,to)mr,p,t),Hm ^ UBLip)] + ^[W{p),SUn{p)], 
iTi — ~ — Ti — — 



(B43) 



(B44) 



where X„ represents the collisions between noncondensate atoms 



P2,P3,P4 

X {l + ffip),i<(p2;P3,P4)} 
-{fflp),A>(P2;P3,P4)} 

— y V- 



fjy2 "^"P ' ^P3 £p4))'^P+P2,P3 + P4 

P2,P3,P4 
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-^P+P2,P3-P4 



P2,P3,P4 
1> 



A (P2;P3,P4) (P2;P3,P4) 



(B46) 



We now comment on the spin rotation terms involving 6U_n and S^,c■ This type of contribution is known as the 
"off-energy-shell" term. In the high-temperature (or low-density) Maxwell-Boltzmann limit T ^ Tc, one can show 
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that 5^_{p) is canceled out by a term second-order in the i-matrix in the mean-field term 0, [o^l ■ In the degenerate 

quantum gas, the role of the ofF-energy-shell term is not obvious, and there are several discussions in literature (see 
discussions in Jeon and MuUin ^18.]). In any event, in our present treatment we do not include renormalization effects 
consistently to order g^. We thus drop the off-energy-shell contribution to the streaming term. For consistency, we 
also drop the anomalous correlation function and the principal-value part of ^ in the GP equation. This means that 

we drop the second-order correction to the condensate chemical potential. Including these second-order contribution 
will require more careful treatment of the many-body perturbation theory. 

Finally, we can replace the momentum sum l/Vj^p by the integral J dp/{2TTfL)^ and the Kronecker delta function 
V^P,p' by the Dirac delta function (27r?i)'^5(p — p'). We then obtain 



xS{ep 



dp2 



dp3 



{2TThf J {2nh) 



dpi 



£pJ)S{p + P2 -P3 -P4) 



{l + ffi(P),A<(P2;P3,P4)} - {mP):A>(P2;P3,P4)} 



(B47) 



IcK(p)] 



dpi 



(27rS)3 



dp2 



dp3 

Pi - P2 



xS{ec + £pi - £p2 - £p3))S{mVc + Pi - P2 - Ps) 
x{Sip ~ pi) [{1 + ffi(pi),4<(C; P2, Pa)} - {^(pi),i>(C; P2, Pa)} 

+(5(p- P2 



-8{p - P3) 



{l + SP2),i<(Pa;C,pi)}-{^(p2),A>(Pa;C,pi)} 
{1 + ffiPa) , i< (P2 ; Pi , C) } - {^(Pa) , i> (P2 ; Pi ; C) } 



(B48) 



R 



dpi 



dp2 



{2TTh)3 J {2TTh) 



dps 



xd{mvc -I- Pi - P2 - Pa)'5(£c + Cpi - £p 
[^^(Pi;P2,Pa) -^^(Pi;P2,Pa)], 



(B49) 



where the various Af- are given in Eq. HB31|I and Eqs. ()ij35|l - HB37|l . 



APPENDIX C: SPIN TRANSPORT RELAXATION TIMES 



In this Appendix, we calculate various spin transport relaxation times that appear in the spin hydrodynamic 
equations. We first consider the following relaxation term: 



dt 



coll 



dpi 



dp2 



dps 



h J {2TThf J {2TThf 
xS{m-Vc + pi - P2 - Ps)S{ec + ip, - ip^ - ip^) 

s,s'=±l 

x{ec + s'e„)(l + 5s,s'){l + fis)f2shs'- 



We then linearize Eq. ICip in ^diff and e„ — Cc to find 



dSn 
dt 



coll 



- I (fin Gc), 



(CI) 



(C2) 



where the two relaxation times are given by 
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xS{pi - P2 - p?,)S{fi.c + Pi - P2 - Pa) 

X^(l+<5,, + )(l + /l,)/2./3+, (C3) 



He 1 Trg'^Uc f dpi f dp^ . 

- 'dp 



n hh J (27rfi)3 J (27rn)3 

xS{pi - P2 - P3)5(/ic + Pl - P2 - Ps) 

X 5^(1 + + /2«)(1 + /3-) - (1 + fu)f2sh-]. (C4) 

S 

We next calculate the relaxation times associated with the spin current. Linearized collision integrals are given by 

SLn = Sl}l + 6Li, (C5) 



where 



SL = Sl}i + Slj, (C6) 



X(5(pi + P2 - P3 - P4)^(£pi + £p2 - ^P3 - £"^4) 

x[(5(p-(5pi) -(5(p-p4)] 

X 1^(1 + • £)(1 + Ss,s'){l + fu){l + f2s')hs'fAs 

s.s' 

a4s + 4s'-4s'-4s), (C7) 



X(5(pi + P2 - P3 - P4)'5(£pi + - £p3 - £p4) 

x[(5(p-pi)-(5(p-p4)] 

X + + f2s')hs'h-s - fufls'O- + hs')0- + /4-.)] 

S.s' 

Xsl- [^i -ifi+ Ss,s'{^2 - V^f ) + Ss,-s'{i^i - ^i)], (C8) 



X(5(pi - P2 - P3)^(Mc +£pi - £p2 - £p3) 

X X](['5(P - Pi) - ^(P - P2) - S{p - P3)]l 

s 

+{s[S{p - Pi) - S{p - P2)] - S{p - P3)}e„ • £) 

X (1 + 5,,+ )(l + /l.)/2«/3+(V'|+ + 4s - 4s)^ (C9) 

x5(pi - p2 - P3)^(a*c + £pi -ep2 - £p3) 

X [(1 + hs)f2-sh+ - hs{l + /2-.)(l + /3+)] 

-<5(P - P3)[(l + hs)hsh- - flsil + /2.)(1 + /3-)] 

X [-v;^ + - vT^) - ■ £). (CIO) 
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When linearizing Xc, we assumed e„ = and ^diff = 0. 

Taking moments of these Unearized coUision integrals with using Eq. H1U5|) , we find 



PTr #111 = T^^pTr [mj, = f^, (CU) 



and 



where the four relaxation times are given by 

1 -ng^ n (3 f dpi f dp2 f dpa 



dp4 



3?i m J {2nh)^ J {2nh)^ J {2nh)^ 

xS{pi + P2 - P3 - P4)'5(£pi + - £p3 - £pj 

X(P1 - P4)'(l + /l + )(l + /2-)/3-/4+, (C13) 

2TTg'^nc h (3 f dpi f dp2 



-j-^ 3h n+h- m J {2iThY J {2iThY 



dp 



3 



X(5(pi - p2 - P3)5(Mc + epi - - ^ps) 
xp2(i + /,_)/2_/3+, (C14) 



1 7rg2 f dpi f dp2 [ dpz , , 

' dpi 



X(5(pi + P2 - P3 - P4)(^(epi + - - ^pa) 

X + + hs')f3s'f4s - fl-sf2s'{l + /3.')(1 + /4s)] 

s,s' 

x[s(pi - Pif + s'(pi - P4) • (pi - P3)], (C15) 

1 ng^ f dpi f dp2 , , 

dp3 



X5(pi - P2 - P3)<5(Aic + Epi - £p2 - Epa) 

X ^{s[(l + h-s)f2s)h+ - h-s{l + /2s)(l + /3+)] 

-[(1 + h^s)f2-sh- - + /2-s)(l + /3-)]} 

4 



x(pi + sp2-P3)}. (C16) 



Finally, the longitudinal spin diffusion time rj^, is given by 

1 1 1 



and the transverse diffusion relaxation time is given by 

1 1 1 



r-L T-L 

^D,n 'D,, 
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